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Abstract 

As is known, tetrahedron equations lead to the commuting family of transfer- 
matrices and provide the integrability of corresponding three-dimensional 
lattice models. We present the modified version of these equations which 
give the commuting family of more complicated two-layer transfer-matrices. 
In the static limit we have succeeded in constructing the solution of these 
equations in terms of elliptic functions. 
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1 Introduction 



Unlike the theory of solvable two-dimensional lattice models of statistical me- 
chanics the theory of three-dimensional models is quite poor. Up to recently 
two-state Zamolodchikov model [||], 0, [Q has been the only nontrivial exam- 
ple of three- dimensional solvable model. In recent paper Bazhanov and 
Baxter have generalized Zamolodchikov model for an arbitrary number of 
states. The Boltzmann weights of their model satisfy the three-dimensional 
star-star equation, tetrahedron equations and possess some remarkable sym- 
metries under the rotations of the elementary cube of the lattice 

i- 

In his first paper [|l| Zamolodchikov have constructed his solution of the 
tetrahedron equations in the so called "static" limit. He interpreted his 
solution in terms of scattering amplitudes of straight strings. The static 
limit corresponds to the scattering of strings with zero velocities. Further 
the static solution was generalized for arbitrary velocities of straight strings 
and all string amplitudes were parameterized by the trigonometric functions 
depending on tetrahedron angles 0. 

However, the first unpublished version of paper [|l|] have contained the 
model in the static limit, whose weight functions were parameterized in terms 
of elliptic functions. It turned out, that these weight functions satisfy only 
some part of tetrahedron equations. The numerical computer tests, carried 
out by Bazhanov in 1980, showed that one can satisfy each concrete tetra- 
hedron equation by changing the signs of weight functions, but all efforts to 
choose the signs of weight functions so that all tetrahedron equations were 
satisfied were futile. 

In this paper we show that the sign factors of weight functions of el- 
liptic Zamolodchikov model can be chosen in a such way that these weight 
functions satisfy the modified system of tetrahedron equations. On obeying 
these equations we have succeeded in constructing the three-dimensional inte- 
grable model with commuting family of two-layer transfer matrices following 
the ideology of paper . 
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2 Weight functions of the model 



The first variant of Zamolodcliikov model of straight strings corresponds to 
the statistical three-dimensional model on a cubic lattice with spin variables 
lying on the faces of the lattice. Baxter (see, for example |^) has shown 
that one can reformulate Zamolodchikov model as an interaction-round-cube 
model with spins in the sites of the lattice. Hereafter we will follow the 
ideology and notations of paper [Q. 

The weight function W{a\efg\bcd\h; 91,62,03) of the model depends on 
the values of eight spin variables, surrounding one elementary cube of the 
lattice (see Fig. 1), and three spectral parameters 61, 6*2, 6*3, satisfying the 
static limit condition: 

61 + 62 + 93 = n. (1) 




= W{a\eJ,g\b,c,d\h;ei,e2^es) 



Fig. 1 



Hereafter all spins take the values 0, 1 and we will imply that all calculations 
over the spin variables are carried out modulo 2. 

The W function satisfies the following "recoloring" symmetry relations: 

W{a\eJ,g\b,c,d\h) = W(a\e, f, g\b,c,d\h) = 
= W{a\e,J,g\b,c,d\h) = W(a\e,J,g\b,c,d\h), (2) 

where we omit the dependence on spectral parameters, a = 1 — a,b = 1 — b, 
etc. Taking into account symmetry relations (^ we have 64 different weight 
functions. 

Define the following combinations of spin variables a, e, f, g, b, c, d, h: 

ji = a + b + e + h, j2 = a + c + f + h, j3 = a + d + g + h, 

nil = e + h, 1712 = f + h, = g + h, (3) 
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and their combinations: 

h = 32 + is, k = ji + is, ^s = ii + 32- (4) 
All spin variables in formulas (H-H) take the values 0, 1 and the weight func- 



tion W in the static limit is defined by formula: 



„J1 „J2 J3 
, 5i *9 5q 



W{a\eJ,g\b,c,d\h;e^,e2,es) = (-lf(n,n,n,rn,,rn,,nr,)j_^3j^nm3^ (5) 

Ci t2 C3 



where we use the following notations: 
sn[ 



1,2,3, (6) 



sn, cn, dn are elliptic functions of modulus k, K is the complete elliptic 
integral of the first kind (see, for example [y]). 

If we take spectral parameters ^j, satisfying the condition < < tt, 
i = 1,2,3, then all values of the elliptic functions in are non-negative, 
and we choose the positive signs of all square roots in formulas (^). 

The sign factor $ in formula (^) is defined by the following expression: 

'^'(ii,i2,i3, mi, 1712, ms) = jim2m3 + j2m3mi + ismim2 + 
+ni{jim2m3 + 32^.37111 + iamimg + 3132^3 + i2i3mi + isiimg) + 
+^2iii2is + J^sO'i + i2 + is) + n^imi + m2 + mg) , (7) 

where ni, n2, n^, are arbitrary integer numbers. 

Hereafter we will suppose static limit condition ([1|) to be valid and use 
the following notation: 

W{a\e, /, g\h, c, d\h; 6^, 62) = W{a\e, f, g\b, c, d\h; 9^, rr - 6^ - 62, 62). (8) 

3 Modified tetrahedron equations 

Let us consider the following equation 

^ 1^(0110102031616263^; 01, 6'2)W(ci|a26362|rfc6C4|64; 6*1, 6*2 + ^s)^ 

d 

xVF(63 1 0402^1616403 1 C5; 6*1 + 6*2, ^3) W^(c^| 6461621030605104; 02, 0-i) = 
= W^(c^|o.2asO'4|o50604|64; 61, 6*2)1^(011^02031610403105; 61, 62 + ^s)x 

d 

Xl^(0i|o2Oi62|0306d|o4; 61 + 02, 6'3)iy(63| 040201 1010203 

M;^2,^3), (9) 
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In appendix we will prove that the weight function W given by formulas 
satisfies modified tetrahedron equations @ with 

W{a\e, /, g\b, c, d\h; 9^, 9^) = W{h\b, c, d\e, f, g\a; 9^, ^2). (10) 

Note that additional weight functions W can be obtained from weight 
function W with the help of the so called {rpY transformation (see 0,[§|]). 
For the Bazhanov-Baxter model, the weight functions are invariant under 
(rp)^ transformation, W = W, and modified tetrahedron equations (||) are 
reduced to the usual tetrahedron equations in the static limit. 

In our case only absolute values of weight functions have this invariance, 
but weights W differ from W by some sign factors. More explicitly, for any 
fixed set of integers rii, n2, ,72,3 ,714 (see formula (|^) we can choose another 
set n'l, n'2, 77,3, 77-4 in a such way that the weight function W is replaced by 
W and vice versa. 

It means that the weight functions W and W satisfy also another equa- 
tion: 

^ W( ai I ci C2 C3 1 61 62 63 ^1 , ^2 ) ( ci I a2 63 1 (ice C4 1 64 ; 6*1 , 6*2 + 6*3 ) X 

d 

xTF(63 1 0402^1616403 1 C5; 9i + 62, 6'3)W^(c^|&46i62|c3C6C5|a4; ^2, 6*3) = 

= X! W {d\a2a^a4\c5CQCi\h4^; 9i, 6'2)W^(ai|(ic2C3 1610403 jcs; 9i, 62 + 6'3)x 
d 

xl^(ci 1020162 |c3C6(i|a4; 9i + 92, 6'3)li^(63|c4C2Ci|oi0203|(i; 92, 93). (11) 



4 Solvable model on three-dimensional lat- 
tice 

As is known, tetrahedron equations lead to the commuting family of transfer- 
matrices and provide the integrability of corresponding three-dimensional 
lattice models [0, [|1^. In paper p it was shown that one can use the 



modified Yang-Baxter equation for constructing the commuting family of 
two-layer transfer-matrices. Similarly, we have succeeded in constructing a 
commuting family of two-layer transfer matrices for the special model on 
three-dimensional cubic lattice, where weight functions W alternate with 
weight functions IV in a checkerboard order in all directions. Note that the 
proof of this fact demands that both equations (pf), (111]) should be valid. 
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Horizontal transfer-matrices T {61,62) constructed with the help of alter- 
nating weights W{6i, 62) and W{6i, 62) (see Fig. 2), satisfy the commutation 
relation: 

[m,^2),m,^3)] = o. (12) 



/ / 




/ / 




/ / 






/ / 
/ / 




/ / 




/ / 



Fig. 2 

More explicitly, we have shown that the weight factor corresponding to the 
composite cube constructed from four "elementary" weights W and four "ele- 
mentary" weights W (see Fig. 3) satisfies the tetrahedron equations of mixed 
type with spin variables placed in the corner sites, in the middles of edges 
and in the centers of all faces of the composite cube. 
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Fig. 3 

The proof of this fact is based on the repeated 16-fold apphcation of modified 
tetrahedron equations (||), ([l^) and will be published in a more detailed 
version of our paper. 

5 Conclusion 

In this paper we have formulated the statistical model on a three-dimensional 
cubic lattice with Boltzmann weights parameterized in terms of elliptic func- 
tions. Further it is naturally to ask about the existence of the deviation 
from the static limit. It seems that the broken (r/9)^-invariance is a principal 
feature of this model. The point is the following: the limit A; ^ of our 
model looks very similar to the static limit of the Zamolodchikov model but 
differs by some sign factors. Hence, the generalization of the static elliptic 
solution (if it exists) cannot be obtained by the deformation of the "full" 
Zamolodchikov solution. 

Also it will be interesting to calculate the statistical sum of this model 
and to investigate its critical behaviour. And at last it is naturally to ask 
about possible a generalization of this model for the number of states N > 2 
by analogy with the Bazhanov-Baxter model. 
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6 Appendix 



In this appendix we give the sketch of the proof of modified tetrahedron 
equations for the weight functions W and W given by formulas d^-H) and 
(P^Ol). First let us introduce some useful notations. Instead of spins a^, bt, 
Ci G Z2 we will use new spin variables with values ±1 and denote them by 
Greek letters: 

ttj = 1 — 2ai, /3i = 1 — 2bi, z = 1, . . . , 4, 

7, = 1-2q, z = 1,...,6. (13) 

Define also nine combinations of spins (|13[): 

Hi = aia2(3iP4, /i2 = aia^PsPi, fi3 = a3a4/32/33, 
Pi = ai/?i7i, P2 = aiPsIs, Ps = "3/^375, P4 = a2/?274, 

ui = q;i/3272, = a4/3i76- (14) 

Also let us introduce the following notations: 

s, = k'/'sn{^-^), c, = ky'p^ ^ = 1,2,3, (15) 



Cl2 = kl . ,n.I,.^u-. ^ C23 = k' - 



5123 = A; /sn( ), ^123 = fc / (16) 

Substituting formulas (||-|D, ([l0|) for W and functions in modified tetra- 
hedron equations (]^) and taking into account formulas (^^-|T^) we obtain 
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(—1)^ 2 2 2 ^ C3 ^ ) 2 X 

. (1-M1M2M3) „^,^ '^^^a-'^ ^l (^l^'ji-^) ,.,. (1 + M3-1--3P1P3) 

1^2 C;^23 C23 J 2 X 

_ (l±AaM2) p2 ^■^i±^^ P2 P3 ~- P2 _il±iilii2M (1-M1P2P3P1P2P3P4) 

l'5l2 '^23 "^2 -5123 j ^ ^ 

(l+Pl) (i±P2l.„ (I+P3) (1+M1P2P3) ^, (1--M2M3) (I-P1M2) 

r-| I ~ 2 Pl ~ 2 P2— 2 P3~ 2 P4— 2 2 w 

[i + /i2Si S2 S3 S123 C12 C23 X 

(1+"2)(1-P3'"1'^3P1P3)(1-A'1M2'"1'^3P2P4) (1+P1P1P2P3P4) | ^^^.^^^^ (1+A'3P1P2P3P4) ] 

_ 1 I ~ 2 Pl~ 2 ^ 1^3PlP3P4~ 3 •^> 3~-- 2 V lP2P3 

— i + /i2Sl S2 S3 S;^23 ^ 

^ (1-M2M3) ^ ^ ^^,^^ (miM2-1) ^ (l+n2)(l-M3^1^3PlP3)(l-MlP2'^l'^3P2P4) 

C12 C23 (-1) 4 X 

^P4KP3 ''+''--^"y2P3P4) +^^^3^3^^ (l+MlPy2P3P4) ]^ 



(17) 



It appears that these relations depend on only nine independent spin 
variables (0). So we have 2^ = 512 separate relations corresponding to the 
different choices of spin variables. We split these relations into the eight 
groups corresponding to the eight possible choices of spins /ij = ±1, z = 
1, 2, 3. Note that some of these equations have the form 

A + B = A + B, A- A = B - B, (18) 

where A and B are the products of elliptic functions. We will call these 
equations as trivial ones. 

First let us consider the case 



/il = 1, /i2 = 1, /i3 = 1. (19) 

It is easily to show that for such choice of fii we have only eight nontrivial 
relations which can be written in the following form: 

~PlP2~P2P3/'1 _|_ ~Pl~P2~P3~PlP2P3\ _ ~Pl ~P3 , ~P2~PlP2P3 (nn\ 

■'12 "'23 l-"- "T ■'I "'2 ■'3 ■'123 J " ''I "'3 "r ■'2 ■'123 • K"^^) 

Detailed analysis shows that for all other choices of spins fii, i = 1,2,3 
equations (0) are reduced to equation (pOl). Note that in some cases we 
should specify the arguments of elliptic functions in equation (^) to some 
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fixed values (for example, 62 + 0^ = vr) and use the quasiperiodic conditions 
for elliptic functions sn(x), cn{x) and dn{x). 

It still remains to check eight relations ( PPD for arbitrary spins pi, p2 and 
P3 taking the values ±1. Note that any choice of spins pi in ( pOD can be 
reduced to the case pi = 1, p2 = I and pa = 1 with the help of the following 
formula: 

sn{iK' + x) = - — , (21) 

k sn[xj 

where K' is the complete integral of the first kind of the complementary 
modulus k' = {1 — k^)2. 

Hence, we should prove equation (|20|) for the only choice pi = 1, p2 = 
1 and ps = 1. This equation can be easily checked by using of addition 
theorems for elliptic functions. It completes our proof of modified tetrahedron 
equations (P). 
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